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In this paper we study the short-time behavior of the Blume-Capel model at the tricritical point 
as well as along the second order critical line. Dynamic and static exponents are estimated by 
exploring scaling relations for the magnetization and its moments at early stage of the dynamic 
, evolution. Our estimates for the dynamic exponent, at the tricritical point, are 2 = 2.215(2) and 
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Numerical simulation in the short-time regime has become an important tool to study phase transitions and critical 
phenomena. The reason is that universality and scaling behavior are already present in the dynamic systems since the 
early stages of their evolution Moreover, this kind of approach reveals the existence of a new and unsuspected 

critical exponent. As shown by Janssen et al. ^ on basis of renormalization group theory, if we tune the parameters 
at their critical values but with initial configurations characterized by nonequilibrium states, the time evolution of 
I ' quantities like magnetization exhibit a polynomial behavior governed by an exponent 9, which is independent of the 
^ , known set of static exponents and of the dynamical critical exponent z. This new exponent characterizes the so 
^sO called "critical initial slip" , the anomalous behavior of the magnetization when the system is quenched to the critical 
. temperature T^. Working with systems without conserved quantities, model A in the terminology of Halperin et al. 
£2 Janssen et al. found a scaling form for the moments of the magnetization, which sets soon after a microscopic 

time scale tmic- Those relations have been confirmed in several numerical experiments |^-^. For the fcth moment of 
the magnetization, this scaling form reads 

§! M'^''\t,T,L,mQ) ^b-''^/''M^'''>{b^''t,b^^''T,b-^L,b='°mo). (1) 

; 

^ , Here b is an arbitrary spatial scaling factor, t is the time evolution and r is the reduced temperature, t ^ {T — Tc)/Tc. 

I ' As usual, the exponents /3 and v are the well-known static exponents, whereas z is the dynamic one. Equation (|l|) 

' ^ ^ depends on the initial magnetization mo and gives origin to the new exponent xq, the scaling dimension of the initial 

g ' magnetization, related to 9 hy 9 — {xq — (3/i')/z. 

^ , For a large lattice size L and small initial magnetization mo, the system in its early stage presents small spatial and 

1| temporal correlation lengths, which may eliminate usual finite size problems. In this limit, if we choose the scaling 
factor b — t^/^ at the critical temperature (r = 0), we obtain 

H' M(i,mo)~moi' (2) 



from the scaling relation (|l|). The exponent 9 has been calculated for the two-dimensional (2d) [g|,[7| and three- 
dimensional (3d) Ising models, 2d 3-state Potts model 1^, Ising model with next-nearest neighbor interactions 
and Ising model with a line of defects | [lO|] . In addition, this short-time universal behavior was found in irreversible 
models with synchronous and continuous time dynamics [Q. In all of those cases, a positive value for 9 has been 
found, that indicates a surprising initial increasing of the magnetization in the short-time regime Iq ~ . This 

effect can be related to a "mean field" behavior since the system presents small correlation length in the beginning of 
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the time evolution. Thus, when the system is quenched to the critical temperature Tc, it feels as being in an ordered 
state since < T*^^ 

On the other hand, as shown by Janssen and Oerding the behavior of a thermodynamic system is more complex 
at a tricritical point (TP) and the corresponding exponent 9 may attain negative values. 

At a tricritical point the magnetization shows a crossover from the logarithmic behavior M(t) ^ mo(ln(t/to)) at 
short times t << TOq^ to t~^^^ power law with logarithmic corrections, M{t) ~ {t/lii{t/to))^^^'^ in 3 dimensions. The 
above behavior can be stated in the generalized form 

M(t) = mo {Ht/to))-''FM (^^) ^ iHt/to))-'' mo j , (3) 

where Fm{x) 1 or Fm{x) ^ 1/x, respectively for vanishing and large arguments. Below 3 dimensions it reduces to 
the scaling form 

Mit) - TOo^^ • (4) 

Here 9 is the exponent related to the tricritical point of the relaxation process at early times and it is expected to 
assume negative values. 

In this paper, we perform short-time Monte Carlo (MC) simulations to explore the critical dynamics of the 2d 
Blume-Capel model. We evaluate the dynamic exponents 9 and z, besides the static exponents v and (3 at the 
tricritical point. To the best of our knowledge this is the first time it is done numerically. We also estimate the 
dynamic exponents along the critical line. We observe a clear trend toward the values of z and 9 for the corresponding 
2d Ising values when the crystal field D becomes large and negative, indicating dynamic universality in the limit 
D -oo. 

The paper is organized as follows. In the next section we present the model and its phase diagram. Sec. Ill contains 
the main scaling relations and describes our short-time MC simulations. Results are presented for critical points on 
the second-order line. In Sec. IV, we explore the short-time dynamics to study the tricritical behavior. Sec. V 
contains a brief outlook and concluding remarks. 



2. THE MODEL 

The Blume-Capel (BC) model is a spin-1 model which has been used to describe the behavior of '^He — ^He 
mixtures along the A line and near the critical mixing point. Apart from its practical interest, the BC model has 
intrinsic interest since it is the simplest generalization of the Ising model (s = 1 /2) exhibiting a rich phase diagram with 
first and second-order transition lines besides a tricritical point. In real systems tricritical points appear in ^He — "^He 
mixtures, such that when a small fraction of ^He is added to ^He, a critical line terminates at a concentration of "^He 
approximately at 0.67. The BC model or its well known generalization, the Blume-Emery-Grifftihs model was 
studied by mean-field approximation, real space and renormalization group schemes Monte Carlo renormalization 
group approach []l9| and finite-size scaling combined with conformal invariance [^-^2| . The Hamiltonian of the two- 
dimensional model is 

H = -J S,Sj+dY,S}. (5) 

<i,j> i=l 

where < i,j > indicates nearest neighbors on lattices and Si = {—1,0,1}. The parameter J is the exchange 
coupling constant and D is the crystal field. We show its phase diagram in Fig. 1. Table I lists points on the 
second order critical line and the tricritical point where we have performed our simulations. Points in Table I were 



obtained from 1 23 1 and from a private communication of the authors in ]22[ | . That table also contains our results for 
the corresponding critical and tricritical exponents. 

We remark that along the critical line, this model presents a critical behavior similar to that of the Ising model. 
However, exactly at the tricritical point the exponents change abruptly. They are given by the dimensions of the 



irreducible representations of the Virasoro algebra ^4 25 1 with central charge (conformal anomaly number) c = 7/10 
pl[ . In [2^ finite-size scaling combined with conformal invariance permited to observe a smooth change between 
Ising-like and tricritical behavior. In finite systems, Ising-like behavior is reached only when D — * —00. In that limit 
77/2(2 — l/i/) — » 0.125 that is the exact value for the Ising model. In our short-time simulations the same kind of 
crossover behavior is observed for the dynamic exponents z and 9 when we move along the critical line. 



Table 1. Critical parameters and exponents for 2d Blume-Capel model. 
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3. NON-EQUILIBRIUM SHORT-TIME DYNAMICS AT A CRITICAL POINT 

In short-time MC simulations critical slowing down can be neglected. It happens because spatial and time correlation 
lengths are small in the early stages of evolution. On the other hand, we need to deal with several samples taken over 
independent initial configurations since the systems which are being simulated are far from equilibrium. In fact this 
approach requires calculation of average (over samples) magnetization and of its moments M'^'^-'(i), 

1 ^= / \ 

^^^'HO = ^E ' (6) 

where (Jij (t) denotes the value of spin i of jth sample at the tth MC sweep. Here denotes the number of samples 
and L'^ is the volume of the system. This kind of simulation is performed Nb times to obtain our final estimates in 
function of t. In this paper, the dynamic evolution of the spins \(Ji\ is local and updated by the heat-bath algorithm. 



A. The critical initial slip 

The evolution of the fcth moment of magnetization in the initial stage of the dynamic relaxation can be obtained from 
Eq. (|l|) for large lattice sizes L at t = with 6 = t^l'- . This yields 

(i, mo) = ^-'^^/''^mC'') (1, i^«/^mo) . (7) 

By expanding the corresponding first moment equation for small mp, we obtain Eq. ^ under the condition that 
mo is sufficiently small, which sets a time scale ~ mg llli^?^ where that phenomena can be observed. 

In Fig. 2a we present our results for the exponent Q at the critical point ksTf,/ J = 1.695 and Dd J — 0, for lattice 
size L = 80 and 5 different initial magnetizations mo. Our estimates for each = 9{mQ) were obtained from Nb — 5 
independent bins with Ns = 10000, for t up to 100 sweeps. Figure 2b illustrates the numerical evaluation of 9 for 
mo = 0.02 from a log-log plot of the magnetization versus time. The linear fitting in Fig. 2a gives d = 0.193(2) with 
goodness of fit ||l Q 0.72. 

Another method has been recently proposed by Tome and de Oliveira to evaluate 9. It avoids the sharp 
preparation of samples with defined and nonzero magnetization and the delicate numerical extrapolation mo — > 0. 
The method is based on the time correlation function of the total magnetization, 



\t=l ]=1 I 



Starting from random initial configurations the above correlation behaves as C{t) ^ , which permits us to obtain 
the exponent 9 from a log-log plot of C{t) versus t. We obtained 9 — 0.194(3) for ksTc/J = 1.695 and Dc/J = 
choosing the time interval [20 — 150] due to the highest value of Q (Q = 0.99). This value is in complete agreement 
with our above estimate of the exponent 9 and it is consistent within error bars with previous results for the 2d Ising 
model. In table I we also present results for 9 at several points of the critical line. 



B. Dynamic critical exponent z 



The observables in short-time analysis are described by different scaling relations according to the initial magnetiza- 
tions. In particular, the second moment M^^^it^V) in Eq. (||), 

\ i=\ I i^j 

with mo = behaves as L"*^ since in the short-time evolution the spatial correlation length is very small when 
compared with the lattice size L. Thus, we arrive at M 



This equation can be used to determine relations involving static critical exponents and the dynamic exponent z 
However, a way to evaluate independently the exponent z is through out the time-dependent fourth-order 
Binder cumulant at the critical temperature (t = 0), 

U4{t,L,ma) = 1 — — -2, (11) 

3(Af(2)(t,L,mo))' 

which obeys the equation 

Ui{t, L, mo) = Ui{h-H, b'^L, 6"^°mo) . (12) 

If we set mo = 0, we eliminate the dependence on the exponent xq and the exponent z can be evaluated through 
scaling collapses of the generalized cumulant for different lattice sizes P,p9|. To match the Binder cumulants Ui{ti, Li) 
and Ui{t2,L2) obtained from two time series for lattice sizes Li and L2, respectively, with b = L2/L1 {L2 > Li), we 
interpolate the series Ui{t,Li) to obtain U4{b~^t, Li). Next, we define the function 

tf 

X\z) = — ^ p4(&-"i, Lr) ~ U4{t, L2)] ' , (13) 
^ ' t=ti 

where the best estimate for z corresponds to the one which minimizes x^(^)- 

In Fig. 3 we show the scaling collapses of the Binder cumulants for different lattice sizes. We have collapsed 
the following pairs of lattices (ii,L2) = (10,20), (20,40) and (40,80). From the largest pair of lattices we obtained 
z = 2.159(6) in the time interval [50 — 1000]. Our final error estimate is based on 25 different collapses obtained from 
Nb — 5 independent bins for each lattice size. 

Another universal behavior of the dynamic relaxation process also described by Eq. (|l|) can be obtained with the 
initial condition mo = 1 ]3C|-|3^ . This condition is related to another fixed point in the context of renormalization 
group approach. Thus, starting from an initial ordered state one obtains a power law decay of the magnetization at 
the critical temperature, 

M{t) - t-f^/"^ , (14) 

when we choose 6~^t = 1 in the limit of L ^ 00. Taking into account this relation, another method has been proposed 
1^ to estimate the dynamic exponent z. This approach uses the second cumulant 

which should take the simple form U2{t, L — > 00) ~ i'*/^. The advantage of this procedure is that curves for different 
lattices lay on the same straight line in a log-log plot without any re-scaling in time. However, this technique has not 
been successful in at least two well known models: the two-dimensional q = 3 Potts model ||^ and the Ising model 
with three spin interactions in just one direction |^ . The reason for the above disagreement may be related to the 
behavior of the second term of r.h.s in Eq. (0) when mo = 1. We have proposed |3^ that this behavior could indeed 
be obtained working with the ratio F2 = M^^^/M^ with different initial conditions for each moment. As we know 



the behavior of the second moment of the magnetization when samples are initiaUy disordered (toq = 0) and also the 
time dependence of the magnetization when samples are initially ordered (mo = 1), we easily obtain 

A log-log plot with error bars for the critical point ksTc/J = 1.695 and Dc/J = is presented in Fig. 4 for L — 160. 
We obtained z = 2.1057(7) with Q = 0.99 in the range [30 — 200], which does not agree with the value obtained from 
Eq. (p!^. However, as we move away from the tricritical point, the values of z obtained (Table 1) with Eq. (|l6|) show 
a clear trend toward the expected value of the dynamic exponent z {z = 2.1567(7)) of the 2d Ising model js^. On the 
other hand, the values obtained from the cumulant in Eq. ( [l^ ) remain essentially the same along the entire critical 
line. It seems the cumulant is less sensitive to crossover effects than the mixed method. 



C. Static exponents and universality class 

The exponent l/vz can be obtained by differentiating lnM(i,T, mg) with respect to the temperature at Tc, 



a In M(i,T, L) 



di 



.l/vz 



(17) 



if we consider the scaling relation for the magnetization when the initial state of samples is ordered (mo = 1) [p9[ . 

Our results for vz were obtained through finite diferences at Tc±6 with 5 = 0.001. They rely on Nb = 5 independent 
bins with Ng = 5000 samples each for L — 160. The time interval [80 — 200] corresponds to the range where the 
goodness of fit parameter attains its highest value {Q = 0.99). 

In Fig. 5 we show the log-log behavior of the derivative 9rlnAf(i) at ksTd J — 1.695 and Dc/J — 0. In Table I 
we present our final estimates of 1/v. The 6th column is obtained with the estimates of z from Eq. ( p^ (data in 4th 
column), while the 7th column corresponds to estimates for v with values of z from Eq. ( p^ (data in 5th column). 

Since we have already collected estimates for vz, it is straightforward to obtain estimates for (3 following Eq. (|lj). 
Our estimates of /3 are presented in the last column. Our values in Table I can be compared with theoretical predictions 
for an Ising like critical point {1/v = l,/3 = 1/8). 



4. RESULTS FROM SHORT-TIME DYNAMICS AT THE TRICRITICAL POINT 



From the results presented in Refs. [jl| and [Q we can describe the time dependence of the magnetization [k = 1) 
for the 2d Blume-Capel model as 



M{t) 



mot'' 



< t < <o , where 9 >0 
to <t <ti 



<t <to , 
to <t <ti 
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critical 
point 

tricritical 
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(18) 



for an initial small magnetization mo. Here ti stands for the time before the system has reached the thermalization. 

We also included in Table I our estimates for 9,z,l/v and (3 at the tricritical point ksTt/J = 0.610 and Dt/J = 
1.9655 working with lattice size L = 80. 

In Fig. 6a we show the values of for 5 different initial magnetizations mo at the tricritical point and in Fig. 6b 
the time evolution of the magnetization for mo — 0.08 in a log-log plot. Our estimates for each 9{mo) were obtained 
from Nb = 20 independent bins with Ns = 10000 samples, for t up to 80 sweeps. The linear extrapolation in Fig. 
6a gives 9 — —0.53(2) with goodness of fit Q = 0.75. The corresponding study with the time correlation function 
in Eq. (||) also gives 9 — —0.53(2) with Q = 0.99 in the time interval [20 — 80]. The generalization of the dynamic 
scaling relation for the ki\i moment of the magnetization at a tricritical point can be written as [p4| 



M(*=) {t, T, g, L, mo) = fe-^ft/^^MW (6"^^, b^/'^r, b'^'/'^g, b-^L, 6^«mo). 



(19) 



It differs from the critical case by the scahng field g that measures the shunting line of critical point along of tangent 
transition line (at the tricritical point). The quantity (pt is the known as the crossover exponent. At the criticality we 
have (j)t — g — 0. 

We show in Fig. 7 scaling collapses of the cumulant U4{t, L) at the tricritical point, whose behavior is quite different 
of the corresponding one at a critical point (Fig. 3). Our first estimate based on Eq. (|l2|) leads to z = 2.21(2). This 
value was obtained from the pair of largest lattice sizes L — 40 and 80. The value remains the same when we consider 
the time interval [10,1000] or [200,1000]. The second estimate for the dynamic exponent, based on Eq. ( [l^ ) gives 
z = 2.215(2). It was obtained from a larger lattice {L = 160) in the time interval [30, 200], with Q = 0.71. We do not 
show the log-log plot of F2{t) in this case because it is quite similar to Fig. 4. We had to restrict the time interval, 
when compared with the Ui calculation, in order to obtain acceptable values for Q. Here, in contrast to the estimates 
of z (4th and 5th column) at the critical points listed in Table I, both methods lead to the same numerical result. 

The evaluation of static exponents v and (3 at the tricritical point follows the same procedure as applied to the 
critical line. Estimates exhibited in Table I are in good agreement with results provided by conformal invariance 
1/ = 9/5 and /3 = 1/24. 

Next, in order to study the influence of the local dynamics over the values of the exponents we recall the simulations 
with Glauber dynamics performed by de Alcantara Bonfim. Our estimates of P/vz obtained from the decay of the 
magnetization (starting from initialy ordered state) is 0.0381(1) with the heat-bath dynamics while (i/i^z = 0.0377(7) 
when updated following Glauber rules. 



5. CONCLUSIONS 



We have performed short-time Monte Carlo simulations to evaluate dynamic and static exponents at critical and 
tricritical points of the spin-1 Blume-Capel model. According to analytical predictions by Janssen and Oerding, a 
negative value for the new exponent 9 was obtained for the tricritical point. The dynamic exponent z was estimated 
by colapsing the fourth-order Binder cumulant and also by the function F2 (t) which explores time evolution of the 
magnetization starting from different initial conditions. The value of z was used to obtain the static exponents (3 and 
v which resulted in good agreement with the exact results provided by the tricritical Ising model [^,Q . The dynamic 
exponents also were calculated along the critical line. Estimates for 9 are in good agreement with known result for the 
2d Ising model. However, near the tricritical point we found strong crossover effects on dynamical exponent z when 
using the recently proposed | p3[ technique based on mixed initial conditions. On the other hand, the fourth-order 
Binder cumulant is almost free from crossover effects. 
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(Fig.l) 1. Phase diagram of the Blume-Capel model. The dashed curve is a first-order transition line and the solid is 
a second order one. These curves are connected by a tricritical point (TP). The marked points (x, •) correspond to the 
simulated values. 



e 



0.192- 



kgT/J= 1.695 




(Fig.2a) 2. Exponent 6 in function of initial magnetizations mo for square lattices with L 
least-square fit to the data. 
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(Fig. 3) 4. Cumulants Ui{t,L) for L = 10,20,40 and 80 for initial magnetization mo = 0. The dots on the lines show 
the cumulants for lattice sizes L/2 rescaled in time with z given by Eq. (|l^). 
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(Fig.4) 5. Time evolution of F2{t) for L — 160 with mixed initial magnetizations. 




(Fig.5) 6. Time evolution of the derivative drinM{t,T)\r=o for L = 160 and initial magnetization mo = 1. 
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(Fig. 6a) 7. Exponent 9 in function of initial magnetizations mo for L = 80 at the tricritical point. The straight line is 
a least-square fit to the data. 
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(Fig.6b) 8. Time evolution of the magnetization for L = 80 and fixed mo — 0.08 at the tricritial point. 
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(Fig.7) 9. Cumulants U4,{t,L) for L = 10,20,40 and 80 for an initial magnetization mo = at the tricritical point. 
The dots on the lines show the cumulants with lattice sizes L/2 rescaled in time with z given by Eq. ([13). 



